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Abstract 

The exponential functional of simple, symmetric random walks with 

negative drift is an infinite polynomial Y = 1 + £1 + £i£2 + £i&£,3 H °f 

independent and identically distributed non-negative random variables. It 
has moments that are rational functions of the variables fik = E(£ fc ) < 1 
with universal coefficients. It turns out that such a coefficient is equal to 
the number of permutations with descent set defined by the multiindex of 
the coefficient. A recursion enumerates all numbers of permutations with 
given descent sets in the form of a Pascal-type triangle. 

1 Introduction 

The present work was induced by the following problem from probability theory, 
cf. [1]. Let (Xj)j?_ 1 be a sequence of independent and identically distributed 

random variables with P(Xj = ±1) = |. Further, let So = 0, S/. = Ylj=i 
(k > 1) be a simple, symmetric random walk. Introduce the following exponen- 
tial functional Y of the random walk with negative drift: 

00 

Y = c MSk -Ai/) = 1 + & + £162 + • ■ • , & = exp(Aj - v), (1) 

fe=0 

where v > 0. Y is an approximation of a widely investigated exponential func- 
tional of Brownian motion, important for example in studying Asian options of 
financial mathematics. To generalize it somewhat, let be independent 

and identically distributed random variables, £j > 0. Consider first the finite 
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polynomials 

= i + £t(i + 6 + 6& + --- + 6---£n) (n>i). 

This equation implies the following equality in distribution (denoted by =): 

Y n = 1 + £y„-i, where £ = £1 and £ is independent of Y n -\- Since F = 

1 + £1 + + + • ■ ■ as n — > 00, in the limit we get a stochastic difference 

equation for the infinite polynomial Y: 

Y = 1 + £Y, 

where £ is independent of Y. 

Then using the binomial theorem, the following recursion can be obtained 
for the pth moment e p = E(Y P ) of Y for any positive integer p: 



1 p- 1 / 
1 IP 



^ \k 



Hk£k, (2) 



supposing fip < 1, where Hk = E(£ fc ), e& — E(Y k ), k > 0. We mention that 
E(y p ) < 00 if and only if \i v < 1, and then /i^ < 1 for any 1 < < p as well, cf. 
[4j Theorem 2]. Observe that ([2]) defines a recursive sequence e p = e p (/ii, . . . , /Lt p ) 
for p > 1 with eo = 1, irrespective of any probability theory background. In 
this recursion the /Zfc's may be considered as variables that may take any value 
except 1 for k > 1. Thus from now on we always suppose that Hk 7^ 1 for k > 1 
and no = 1. 

It follows from @ by induction that e p is a rational function of the variables 
Hi, ■ . ■ , Hp f° r anv integer p > 1: 

^ (i-.o-U-m E 4 P) .,.X ( 3 ) 

11 Wj 11 Mpj Wi,..A-0€{o,i}--i 

where the coefficients of the numerator are universal constants, being indepen- 
dent of the parameters in the recursion. 

(v) 

These coefficients j make a symmetrical, Pascal's triangle- like table 
if each row is listed in the increasing order of the binary numbers j p _i2 p ~ 2 + 
• • • + ji2°, defined by the multiindices (ji, . . . , see the rows p — 1, . . . , 5 

in Table 1. 

Two natural questions may arise at this point, independently of any prob- 
ability theory background mentioned above. First, suppose that one defines a 
recursive sequence (e p )^L 1 by ^) with coefficients ■ _ l given by ©. Can 

one attach any direct mathematical meaning to these coefficients „■ then? 
The answer is yes, and rather surprisingly (as was conjectured in [J]), the coef- 
ficient a/j^ j i is equal to the number of permutations 7r £ S p having descent 
7r(i) > 7r(i + 1) exactly where ji = 1, 1 < i < p — 1, cf. Theorem [TJ below. 

Second, can one give a direct way to evaluate the coefficients of „ ? The 

. . . . . Ox—Jp-i 

affirmative answer to this question is partly included in the previous answer, 
since several formulae are known for the number of permutations with given 
descent sets. However, an apparently new recursion was conjectured in j4] , which 
is analogous to the recursion of binomial coefficients in the ordinary Pascal's 
triangle. The proof of this is the content of Lemma [5] below. 
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Table 1: The Pascal's triangle-like table of the coefficients 
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2 The results 

In the next lemma we establish a direct recursion for the coefficients a?f ■ 

3\—3p-l 

Lemma 1. Fix a multiindex (Ji, . . . , j p -i) £ {0,1} P_1 . Let S be the set of 

indices k where jk = 1: 

p-l 

S = {si,... ,s m } = {k : jk = 1, 1 < k <p- 1}, = (4) 

fe=i 

Then the coefficient j i defined by ^j) can be obtained by the recursion 



, - e(^>< '" ,A ' 

fc=0 



1=0 



E if (-ir-'oj^^, (5 ) 



where, by definition, = 1, jo = 1, Sq = and —1 powered to an empty sum 
is 1. 

Proof. The second equality in ([5]) is a direct consequence of the definitions above. 
To show the first equality, substitute into ©: 

6P ~ l-[l p ^ Ui(l- Ml )---(1- M ) ^ Sl.^-^1 Mfc-l- 

Here, multiplying by the common denominator and then collecting the coeffi- 
cients of fJ. 3 i ■ ■ ■ MpTi for each (ji, . . . , we obtain 



s P (l - Ml) •••(! - Mp) 

.A- 



fc =o 0'i,...,j fc -i)e{o,i}'=- 



p-i 



E >'- ; •••/';; ,E I ^,-.^-^ 1 • • • c-^- 1 - 

(ii,...,i P -i)e{o,i}p-! fc=o v J 

This and ([3]) imply the first equality in ([3]). □ 
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Now we turn to the proof of the equality of the coefficient ^ _ t given by 
© and the number of permutations b^' (S) with descent set S given by (@| . The 
descent set of a permutation tt € S p is defined as D(tt) — {i : 7r(i) > 7r(i+l), 1 < 
i < p — 1}. It is known, cf. [3j p. 69], that the number of permutations it *E S p 
with a given descent set S = (s%, . . . , s m ), 1 < si < • • • < s m < p — 1, can be 
obtained by the following inclusion-exclusion formula: 

b ( P )( S) = 6 (p)( Sl ,... ;Sm ) 



l<ii < ■ ■ -<ij <m 



Theorem 1. TTie coefficient j _ 1 given by HP and (0) is eguaZ to £/ie number 
of permutations b^ p \S) with descent set S given by Q). 

Proof. It is enough to show that the numbers b^ p ' (S) satisfy the same recursion 

(T3T) as the numbers oi p ' • do, that is, 
* ' ji-'-jp— l 

fcW( Sl , . . . , s m ) = ]T r (-1)™- ; 6^)( S1 , . . . , *,_!), (7) 

where, by definition, so = and 6( s °- ) = fr^ - 1 = 1. 

To show this, let us substitute the inclusion-exclusion formula ((SJ) into the 
right hand side of ([7]): 

rn / \ 

= H) m +EL Hr' E 

;=1 \ 1<£l<---<V<2-1 
rn 

= (-i) m + E E (-i)™-*" 1 

1=1 l<i 1 <-<i j <l 



1-1-3 



Si 

'«! ) ^2 ? * * * 1 "J ^~ 



21 ? ^1 ) ' * ' ) Si - , J* 

m-(j + l) / P 



l<«i <---<i 3 <l<m J 



= (-ir+ E 

l<2i<---<ij<Z<m 

E (-ir~ r ( p _ ) =b^( Sl ,...,, 

l<t 1 <---<i r < m \ S h> S i2 s n ,...,s tr 8i r _ lt p s lr j 

This proves ([7]), and so completes the proof. □ 

Lemma [JJ described a recursion that uses all previous rows of Table 1 to 
compose the elements of a new row. In the next lemma we show a recursion that 
uses only the previous row and which is an analog of the recursion formula in the 
ordinary Pascal's triangle: = + C^ 1 ) ■ Interestingly, an application of 

this identity is a key step in the following algebraic proof as well. We also give 
a simple combinatorial proof which basically translates the well-known method 
by which permutations in S p can be obtained from permutations in S p -i by 
adjoining the number p. 
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Lemma 2. The following recursion holds for anyp > 2 and multiindex (ji, . . . ,j P -i) E 
{0,1} p 



,(p) 



5> 



a 




(p-i) 



E 



« 



(p-i) 

'il...ip_ 2 ' 



(8) 



(i 1 ,...,»p_ 2 )6i(ji)— Jp-i) 



w/iere = 1, 5, = \j t - j»_x | /or i > 2, <S X = 1, = jk for 1 < k < i - 1, 
= jfe+i /or i < k < p — 2, and L(ji, . . . ,j P -i) is the set of all distinct 
binary sequences obtained from (ji, . . . , j p -i) by deleting exactly one digit. For 

example, a$ 10 = 11 = a$ + a$ + a$ x . 

Proof. First we prove the second equality in ([8]). For this it is enough to show 
that if the same binary sequence is obtained from (ji, . . . ,j p —i) when eliminating 
either the fcth or the Zth digit (k < I), then all digits between the fcth and 
Zth (including these two) are uniformly either O's or l's (a run of O's or l's). 
Therefore, the two recursions given in ([5]) are the same. 

Consider a multiindex (ji, . . . ,jk-i, jk, ■ ■ ■ ji+i, ■ ■ ■ , jp-i) E {0, lj^^ 1 . 
Suppose that we get the same binary sequence by deleting jk and ji, respectively: 

(ji, • ■ ■ ,3k-i,jk+i, ■ ■ ■ ■ ■ ,jp-i) = (ji,- ■ ■ ,jk-i,jk, ■ ■ -,3l-i,jl+i, ■ ■ -,jp-i)- 

Then jk = jk+i = ■ ■ ■ = = ji, so the second equality in §5§ really holds. 

Now it remains to show the first equality in ([8|). that is, the recursion itself. 

A combinatorial proof of the recursion. Given a binary sequence (ji, . . . , j p —i ) , 
let us remove a single 1 from a run of l's or a single from a run of O's. Count 
the number of permutations in S p -\ determined by the resulting multiindex 

(ii,... , ip-i). This number is J_ by Theorem [TJ We want to show that 
there is a uniquely determined adjoining of the number p to any such permuta- 
tion from Sp-i to obtain a permutation from S p corresponding to the original 
multindex (j u j p _i). 

If a was deleted from a run of O's, the number p should be inserted right 
after the number at the position of the first 1 following the affected run of O's. 
(If the given run happens to be the last, p is inserted as the last number.) When 
a 1 was deleted from a run of l's, the number p should be inserted right after 
the number at the position of the last preceding the affected run of l's. (If 
the given run happens to be the first, p is inserted as the first number.) Since 
these insertions are the only ones that reconstruct the original descent set, the 
recursion is proved. 

An algebraic proof of the recursion. We are going to proceed by induction 
over p. Thus suppose that the recursion holds for all multiindices of lengths 
smaller than p — 1. First we use the recursion of Lemma [T] for the terms in the 
right side of ([%])■ then we change the order of summation to obtain 




a 



i=l 



■5 



Here in the last expression, one may use the induction hypothesis for the 
first sum. In the second sum observe that l)- 5 * is I if = 1 and ji = 0, 
it is —I if ji—i = and ji = 1, and it equals otherwise. Hence we get the 
identity jk Y^=k+i 1) J< = ifc(l — ip— l)- Thus one obtains that 

X>^« =E{( p ; 1 )i^(-iy- + ™ + *- i -S31 

i=i Jl - 3p - 2 fc =o fc / 

+( p fc 1 )(l-j p -l)j fe (-l) 

- E{(rOK p ;>-^4* ( - 1 ^^ 10 ^ 

To rewrite the terms above we use recursion ((SJ in the following case: 

plus the identity — lpp- 1 = j p -i, and the recursion for binomial coeffi- 

cients: 

+d - j P -i)(-iy i+ -^- + (J 1 2)^-1^:!, 
= E^^c- 1 )^ 1 ^- 1 ^-^^..^- 

This completes the proof. □ 

The results above imply that Table 1 has properties analogous to the ones 

(p) 

of Pascal's triangle: each entry a)?' ^ is a positive integer, the first and the 

last entries, a^o an< ^ a i'\ are 1) the table has symmetries a^' ^ _ x = 

and „■ = of \ -, „■ , and the sum of the 2 P ~ 1 entries in the pth row is 
pi. 

3 Some remarks 

We mention that beside the indirect recursion ([2]) , ([3]) , direct recursions ([5]) and 
©, and sieve formula ©, there are other methods as well in the literature 
that can be used for evaluating the coefficients ■ . Here we mention two 
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of them. Zabrocki [5] uses the following rather fast and practical 'splitting' 
recursion: 



n (p) (p-i\tk-i) jp-k) 

u h---ip-i l^i \j. | j /: .• , r 

{fc:jh-l=0,j*=l,l<*<p} V 7 

where, by definition, jo = and j p = 1. Its verification is simple: this recursion 
divides the set of permutations in S p with descent set defined by the l's in 
(ji, . . . ,j p —i) into disjoint subsets, where the largest number p is at the different 
local maxima k where jk-i = and jk = 1< Then there are (?Zi) ways to pick 

the numbers to precede the largest, with a[ k permutations, and af k ^ , 

r . Ji"-Jk-a ' Jfe+i---Jp-i 

permutations for the numbers succeeding the largest, fitting the given descent 
set. 

Another method (of mostly theoretical interest) can be based on MacMa- 
hon's determinant, cf. [31 p. 69]: 



a 



(p) 

ji—jp-i 



6 (p) (si,...,s m ) = p\ det [l/(s j+1 -Sj)!], G [0,m] x [0,m]., 



where the descent set S = (si, . . . , s m ) is defined by ([!]), so = 0, and s m +i = p. 
This gives the recursion 

6 (p) (si, . . . ,s m ) = -^6 (p) (s 2 - si, . . . ,s m - si) ^r6 (p) (s3 - s 2 , • . . , s m - s 2 ). 

si! s 2 ! 

Obviously, the coefficients a)? j are closely related to other important 
classifications of permutations as well. The Eulerian number A(p, k) which 
counts the permutations in S p having exactly k — 1 descents (that is, exactly k 
runs), where p > 1,1 < k < p, can be written as 

{(ju---jp-i)£{0,l}P-i-.-£ P iZln=l<:-l} 

Also, let I(p, k) be the number of permutations in S p with k inversions 
(p > 1,0 < fc < (I)). By MacMahon's theorem, cf. [U Section 5.1.1], I(p,k) is 
equal to the number of permutations in S p with major index k, thus 



ji—jp-i ' 
{(ji,--^ P -i)e{0,l}*-i : £*r 1 1 i»=fc} 



Finally, let us express the generating function of the coefficients cSj^ ^ _ ± by 
the help of a suitable recursive sequence {e p ) p a =1 given by @. First, one can 
assign a positive integer n to any p > 1 and multiindex (ji, . . . , j p -i) by the 
equation n = 2 P_1 + 53fe=i ifc2 fc_1 , in a one-to-one way. Then let us introduce 



the notation a„ = ■ _ 1 , n > 1. This way the coefficients are arranged in a 
single sequence. 

Second, take the variables (moments) fXk — x 2 for k > 1, x ^ 1, and let 
= 1- Then define the sequence (e p (x)) p *L 1 by ([2]): 
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Also, using © and the definition of n one obtains that 



1 



2 P -1 



e p (x) 



(1 -x)(l -x 2 ){l - a; 4 )---(l -x- 



} J a n x 



n 



Hence the generating function of the sequence (a„)™ =1 can be expressed as 



References 

[1] Knuth, D.E. (1973) The Art of Computer Programming. Vol. 3. Sorting and 
searching. Addison- Wesley, Reading, Mass. 

[2] MacMahon, P.A. (1915) Combinatory Analysis. Vol. 1. Cambridge Univ. 
Press, Cambridge, England. 

[3] Stanley, R. (1986) Enumerative Combinatorics. Vol. 1. Wadsworth and 
Brooks/Cole Mathematics Series, Monterey, Calif. 

[4] Szabados, T. and Szekely, B. (2003) An exponential functional of random 
walks. J. Appl. Prob. to appear. 

[5] Zabrocki, M. (2001) Integer sequence A060351. On-line encyclopedia of inte- 
ger sequences. |http://www.research.att.com/~njas/seque nces/Seis.html 



oo 



oc 




£ x 2 " \l - x)(l - x 2 )(l - x A ) ■ ■ ■ (1 - x 2 " l )e p {x). 
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